How large should whales be?

The evolution and distribution of species body sizes for terrestrial mammals is well-explained by a macroevolutionary tradeoff between short-term selective advantages and long-term extinction risks from increased species body size, unfolding above the 2g minimum size induced by thermoregulation in air. Here, we consider whether this same tradeoff, formalized as a constrained convection-reaction-diffusion system, can also explain the sizes of fully aquatic mammals, which have not previously been considered. By replacing the terrestrial minimum with a pelagic one, at roughly 7000g, the terrestrial mammal tradeoff model accurately predicts, with no tunable parameters, the observed body masses of all extant cetacean species, including the 175,000,000g Blue Whale. This strong agreement between theory and data suggests that a universal macroevolutionary tradeoff governs body size evolution for all mammals, regardless of their habitat. The dramatic sizes of cetaceans can thus be attributed mainly to the increased convective heat loss is water, which shifts the species size distribution upward and pushes its right tail into ranges inaccessible to terrestrial mammals. Under this macroevolutionary tradeoff, the largest expected species occurs where the rate at which smaller-bodied species move up into large-bodied niches approximately equals the rate at which extinction removes them.Comment: 7 pages, 3 figures, 2 data table

Bayesian approaches to reverse engineer cellular systems: a simulation study on nonlinear Gaussian networks

BACKGROUND. Reverse engineering cellular networks is currently one of the most challenging problems in systems biology. Dynamic Bayesian networks (DBNs) seem to be particularly suitable for inferring relationships between cellular variables from the analysis of time series measurements of mRNA or protein concentrations. As evaluating inference results on a real dataset is controversial, the use of simulated data has been proposed. However, DBN approaches that use continuous variables, thus avoiding the information loss associated with discretization, have not yet been extensively assessed, and most of the proposed approaches have dealt with linear Gaussian models. RESULTS. We propose a generalization of dynamic Gaussian networks to accommodate nonlinear dependencies between variables. As a benchmark dataset to test the new approach, we used data from a mathematical model of cell cycle control in budding yeast that realistically reproduces the complexity of a cellular system. We evaluated the ability of the networks to describe the dynamics of cellular systems and their precision in reconstructing the true underlying causal relationships between variables. We also tested the robustness of the results by analyzing the effect of noise on the data, and the impact of a different sampling time. CONCLUSION. The results confirmed that DBNs with Gaussian models can be effectively exploited for a first level analysis of data from complex cellular systems. The inferred models are parsimonious and have a satisfying goodness of fit. Furthermore, the networks not only offer a phenomenological description of the dynamics of cellular systems, but are also able to suggest hypotheses concerning the causal interactions between variables. The proposed nonlinear generalization of Gaussian models yielded models characterized by a slightly lower goodness of fit than the linear model, but a better ability to recover the true underlying connections between variables.Italian Ministry of University and Scientific Research; National Institutes of Health & National Human Genome Research Institute (HG003354-01A2); Collegio Ghislieri, Pavia Italy fellowshi

A Relative Variation-Based Method to Unraveling Gene Regulatory Networks

Gene regulatory network (GRN) reconstruction is essential in understanding the functioning and pathology of a biological system. Extensive models and algorithms have been developed to unravel a GRN. The DREAM project aims to clarify both advantages and disadvantages of these methods from an application viewpoint. An interesting yet surprising observation is that compared with complicated methods like those based on nonlinear differential equations, etc., methods based on a simple statistics, such as the so-called -score, usually perform better. A fundamental problem with the -score, however, is that direct and indirect regulations can not be easily distinguished. To overcome this drawback, a relative expression level variation (RELV) based GRN inference algorithm is suggested in this paper, which consists of three major steps. Firstly, on the basis of wild type and single gene knockout/knockdown experimental data, the magnitude of RELV of a gene is estimated. Secondly, probability for the existence of a direct regulation from a perturbed gene to a measured gene is estimated, which is further utilized to estimate whether a gene can be regulated by other genes. Finally, the normalized RELVs are modified to make genes with an estimated zero in-degree have smaller RELVs in magnitude than the other genes, which is used afterwards in queuing possibilities of the existence of direct regulations among genes and therefore leads to an estimate on the GRN topology. This method can in principle avoid the so-called cascade errors under certain situations. Computational results with the Size 100 sub-challenges of DREAM3 and DREAM4 show that, compared with the -score based method, prediction performances can be substantially improved, especially the AUPR specification. Moreover, it can even outperform the best team of both DREAM3 and DREAM4. Furthermore, the high precision of the obtained most reliable predictions shows that the suggested algorithm may be very helpful in guiding biological experiment designs

Determination of glucose exchange rates and permeability of erythrocyte membrane in preeclampsia and subsequent oxidative stress-related protein damage using dynamic-19F-NMR

The cause of the pregnancy condition preeclampsia (PE) is thought to be endothelial dysfunction caused by oxidative stress. As abnormal glucose tolerance has also been associated with PE, we use a fluorinated-mimic of this metabolite to establish whether any oxidative damage to lipids and proteins in the erythrocyte membrane has increased cell membrane permeability. Data were acquired using 19F Dynamic-NMR (DNMR) to measure exchange of 3-fluoro-3-deoxyglucose (3-FDG) across the membrane of erythrocytes from 10 pregnant women (5 healthy control women, and 5 from women suffering from PE). Magnetisation transfer was measured using the 1D selective inversion and 2D EXSY pulse sequences, over a range of time delays. Integrated intensities from these experiments were used in matrix diagonalisation to estimate the values of the rate constants of exchange and membrane permeability. No significant differences were observed for the rate of exchange of 3-FDG and membrane permeability between healthy pregnant women and those suffering from PE, leading us to conclude that no oxidative damage had occurred at this carrier-protein site in the membrane

Shared component modelling as an alternative to assess geographical variations in medical practice: gender inequalities in hospital admissions for chronic diseases

<p>Abstract</p> <p>Background</p> <p>Small area analysis is the most prevalent methodological approach in the study of unwarranted and systematic variation in medical practice at geographical level. Several of its limitations drive researchers to use disease mapping methods -deemed as a valuable alternative. This work aims at exploring these techniques using - as a case of study- the gender differences in rates of hospitalization in elderly patients with chronic diseases.</p> <p>Methods</p> <p>Design and study setting: An empirical study of 538,358 hospitalizations affecting individuals aged over 75, who were admitted due to a chronic condition in 2006, were used to compare Small Area Analysis (SAVA), the Besag-York-Mollie (BYM) modelling and the Shared Component Modelling (SCM). Main endpoint: Gender spatial variation was measured, as follows: SAVA estimated gender-specific utilization ratio; BYM estimated the fraction of variance attributable to spatial correlation in each gender; and, SCM estimated the fraction of variance shared by the two genders, and those specific for each one.</p> <p>Results</p> <p>Hospitalization rates due to chronic diseases in the elderly were higher in men (median per area 21.4 per 100 inhabitants, interquartile range: 17.6 to 25.0) than in women (median per area 13.7 per 100, interquartile range: 10.8 to 16.6). Whereas Utilization Ratios showed a similar geographical pattern of variation in both genders, BYM found a high fraction of variation attributable to spatial correlation in both men (71%, CI95%: 50 to 94) and women (62%, CI95%: 45 to 77). In turn, SCM showed that the geographical admission pattern was mainly shared, with just 6% (CI95%: 4 to 8) of variation specific to the women component.</p> <p>Conclusions</p> <p>Whereas SAVA and BYM focused on the magnitude of variation and on allocating where variability cannot be due to chance, SCM signalled discrepant areas where latent factors would differently affect men and women.</p

The economic pressures for biosimilar drug use in cancer medicine

The main rationale for using biosimilar drugs is for cost saving. The market development for biosimilar drugs will therefore depend on the degree to which cost saving measures are required by nations, medical insurers and individuals and the absolute savings that could be gained by switching from original drugs. This paper is designed to discover the degree to which financial constraints will drive future health spending and to discover if legal or safety issues could impact on any trend. A structured literature search was performed for papers and documents to 27 August 2011. Where multiple sources of data were available on a topic, data from papers and reports by multinational or national bodies were used in preference to data from regions or individual hospitals. Almost all health systems face current significant cost pressures. The twin driver of increasing cancer prevalence as populations age and cancer medicine costs rising faster than inflation places oncology as the most significant single cost problem. For some countries, this is predicted to make medicine unaffordable within a decade. Most developed countries have planned to embrace biosimilar use as a cost-control measure. Biosimilar introduction into the EU has already forced prices down, both the price of biosimilar drugs and competitive price reductions in originator drugs. Compound annual growth rates of use have been predicted at 65.8% per year. Most developed countries have planned to embrace biosimilar use as a major cost-control measure. Only legal blocks and safety concerns are likely to act against this trend. For centralised healthcare systems, and those with a strong tradition of generic medicine use, biosimilar use will clearly rise with predictions of more than 80% of prescriptions of some biologic drugs within 1 year of market entry in the USA. Delaying the implementation of such programmes however risks a real crisis in healthcare delivery for many countries and hospitals that few can now afford

Inferring transcriptional compensation interactions in yeast via stepwise structure equation modeling

<p>Abstract</p> <p>Background</p> <p>With the abundant information produced by microarray technology, various approaches have been proposed to infer transcriptional regulatory networks. However, few approaches have studied subtle and indirect interaction such as genetic compensation, the existence of which is widely recognized although its mechanism has yet to be clarified. Furthermore, when inferring gene networks most models include only observed variables whereas latent factors, such as proteins and mRNA degradation that are not measured by microarrays, do participate in networks in reality.</p> <p>Results</p> <p>Motivated by inferring transcriptional compensation (TC) interactions in yeast, a stepwise structural equation modeling algorithm (SSEM) is developed. In addition to observed variables, SSEM also incorporates hidden variables to capture interactions (or regulations) from latent factors. Simulated gene networks are used to determine with which of six possible model selection criteria (MSC) SSEM works best. SSEM with Bayesian information criterion (BIC) results in the highest true positive rates, the largest percentage of correctly predicted interactions from all existing interactions, and the highest true negative (non-existing interactions) rates. Next, we apply SSEM using real microarray data to infer TC interactions among (1) small groups of genes that are synthetic sick or lethal (SSL) to SGS1, and (2) a group of SSL pairs of 51 yeast genes involved in DNA synthesis and repair that are of interest. For (1), SSEM with BIC is shown to outperform three Bayesian network algorithms and a multivariate autoregressive model, checked against the results of qRT-PCR experiments. The predictions for (2) are shown to coincide with several known pathways of Sgs1 and its partners that are involved in DNA replication, recombination and repair. In addition, experimentally testable interactions of Rad27 are predicted.</p> <p>Conclusion</p> <p>SSEM is a useful tool for inferring genetic networks, and the results reinforce the possibility of predicting pathways of protein complexes via genetic interactions.</p

Characterizing Dynamic Changes in the Human Blood Transcriptional Network

Gene expression data generated systematically in a given system over multiple time points provides a source of perturbation that can be leveraged to infer causal relationships among genes explaining network changes. Previously, we showed that food intake has a large impact on blood gene expression patterns and that these responses, either in terms of gene expression level or gene-gene connectivity, are strongly associated with metabolic diseases. In this study, we explored which genes drive the changes of gene expression patterns in response to time and food intake. We applied the Granger causality test and the dynamic Bayesian network to gene expression data generated from blood samples collected at multiple time points during the course of a day. The simulation result shows that combining many short time series together is as powerful to infer Granger causality as using a single long time series. Using the Granger causality test, we identified genes that were supported as the most likely causal candidates for the coordinated temporal changes in the network. These results show that PER1 is a key regulator of the blood transcriptional network, in which multiple biological processes are under circadian rhythm regulation. The fasted and fed dynamic Bayesian networks showed that over 72% of dynamic connections are self links. Finally, we show that different processes such as inflammation and lipid metabolism, which are disconnected in the static network, become dynamically linked in response to food intake, which would suggest that increasing nutritional load leads to coordinate regulation of these biological processes. In conclusion, our results suggest that food intake has a profound impact on the dynamic co-regulation of multiple biological processes, such as metabolism, immune response, apoptosis and circadian rhythm. The results could have broader implications for the design of studies of disease association and drug response in clinical trials

Casual Compressive Sensing for Gene Network Inference

We propose a novel framework for studying causal inference of gene interactions using a combination of compressive sensing and Granger causality techniques. The gist of the approach is to discover sparse linear dependencies between time series of gene expressions via a Granger-type elimination method. The method is tested on the Gardner dataset for the SOS network in E. coli, for which both known and unknown causal relationships are discovered